TPI-MINN-98/06-T 
UMN-TH- 1702/98 
ICN-UNAM 98-04 
May 1998 
hep-th/98UoTO 



Energy Reflection Symmetry of Lie- Algebraic 
Problems: Where the Quasiclassical and Weak 
Coupling Expansions Meet 



M. Shifman 

Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 54555 

and 
A. Turbiner 

Institute de Ciencias Nucleares, UNAM, 04510 Mexico D.F., Mexico 

and 

ITEP, Moscow 117259, Russia f 
Abstract 

We construct a class of one-dimensional Lie-algebraic problems based on sl(2) 
where the spectrum in the algebraic sector has a dynamical symmetry E <-> —E. 
All 2j + 1 eigenfunctions in the algebraic sector are paired, and inside each pair 
are related to each other by a simple analytic continuation x — > ix, except the zero 
mode appearing if j is integer. At j ' ^ oo the energy of the highest level in the 
algebraic sector can be calculated by virtue of the quasiclassical expansion, while 
the energy of the ground state can be calculated as a weak coupling expansion. The 
both series coincide identically. 



On leave of absence. 



1. Introduction. A hidden algebraic structure of the quasi-exactly solvable 
(QES) Hamiltonians [U |4j] leads to non-trivial dynamical properties of the QES 
systems. One of such properties was observed in 0: it was noted that all levels in 
the algebraic sector of the simplest QES problem (see Eq. ([!]) below) are symmetric 
under E <-> —E. This property will be referred to as the energy-reflection (ER) sym- 
metry. In the present paper we derive a class of one-dimensional QES Hamiltonian 
with this property and explore the consequences of the ER symmetry. A relation 
between the weak coupling and quasiclassical expansions will be established. 

One-dimensional QES problems are based on a (hidden) si (2) algebra; they are 
characterized by one quantized (cohomology) parameter j where j is half-integer 
||. [7J. The number of levels in the algebraic sector is 2j + 1. In the systems to 
be constructed below each state in the algebraic sector with the energy eigenvalue 
—E (E > 0) is accompanied by a counterpart with energy E, if 2j + 1 is even. 
If 2j + 1 is odd, a zero mode exists while the remaining 2j levels come in pairs 
{ip-E, iPe}- The eigenfunctions of the ER-symmetric levels are related to each other 
by a straightforward analytic continuation 

At large j the number of states in the algebraic sector is large. The highest 
levels still belonging to the algebraic sector can be regarded as highly excited states, 
and as such, are amenable to the quasiclassical treatment ||. The parameter of the 
quasiclassical expansion is At the same time, the lowest levels from the algebraic 
sector are close to those of the harmonic oscillator. The anharmonicity is small, and 
is determined by a small parameter related to 1/j. Under the circumstances, one can 
develop a standard weak-coupling perturbation theory and calculate E as a series 
in the weak coupling. Since the energy eigenvalues of the highly excited and low- 
lying ER-partners coincide, up to sign, the quasiclassical expansion and the weak 
coupling expansion in the QES problems with the ER symmetry must be identical. 
We discuss how this identity is implemented, taking as a representative example the 
ground state and its counterpart. 

The simplest QES problem with the ER symmetry known for a long time || is 
the sextic anharmonic oscillator, with a quantized coefficient in front of x 2 , 
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p 2 + (x 6 -(8j + 3)x 2 )} , j = 0,i 1,~ ... (1) 



In this case the algebraic sector consists of 2j+l levels of positive parity. The general 
QES potential possessing the ER symmetry involves certain elliptic functions and 
is related to some problems of practical importance. 

2. Generalities. The strategy we follow is described in n (see also [|7|]) while 
the notation is borrowed from ||. The generators of the sl(2) algebra are defined 
as follows: 

T + = 2jt;-e^, T° = -j + ^ T~ = ±. (2) 
d£ d£ d£ 
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If j is a non-negative half-integer number, a finite-dimensional irreducible represen- 
tation exists, 

= {?,?, (3) 

where the subscript indicates the dimension of the representation. In general, the 
generators have the meaning of the raising (lowering) operators, 

Rn — > R n ±i , while R n — > R n . 

The generic QES sZ(2)-based Hamiltonian is representable as a quadratic combina- 
tion of the generators T^ 1 and T°, 

H = (C ab T a T b + C a T a ) +C, (4) 
±,o 

where C a b, C a , C are parameters. One can always get rid of C + and C_ in favor of 
C +0 and C _, respectively, due to the si (2) commutation relations. Moreover, C + _ 

and C (. can be eliminated in favor of Coo, as a consequence of the irreducibility of 

the representation Rzj+i (i-e. T + T _ + T~T + + 2T°T° = 2j(j + 1)). The reference 
point for the energy is fixed by putting (7=0. 

After a change of variable and a (quasi)gauge transformation the operator H can 
be always reduced to the Schrodinger form 

H^H = e-^He^\ =ax) = + V(x) . (5) 

The key element in constructing the QES Hamiltonians with the ER symmetry is 
the following observation |§: any tridiagonal matrix of the form 
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(6) 



leads to the characteristic equation 

EP n/2 (E 2 ) = , (n even), and P (n+1 y 2 (E 2 ) = , (n odd) , (7) 

where P n / 2 (z) and P( n +i)/2 (z) are polynomials of z of degree n/2 and (n + l)/2, 
respectively, and n is defined in Eq. (|G]). Thus, any matrix of the form ([]) guarantees 
the ER symmetry of the spectrum. It is evident, that the Lie- algebraic Hamiltonian 
(||) has the matrix representation (|6]) provided that the sum in Eq. (^) does not 
include the terms T + T + , T _ T _ , T°T° and T°. Thus, the most general form of 
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the Lie-algebraic Hamiltonian compatible with Eq. (0) - which, as was explained, 
ensures the ER symmetry in the algebraic sector - is Q 



H = aT + T° + (3T°T- +-/T + + 5T- = 

m£ + B®± + C®, (8) 

where a, /?, 7 and 5 are numerical constants, and A, B, C are polynomials in £ of the 
third, second and first degree, respectively, 

A(£) = -af+ft, B(() = [a(3j-l)-j]e + (5-Pj) 7 C(0 = 2( 7 j -af)£ . (9) 

Not all of the four constants above represent physically interesting parameters. In 
general, two constants can be fixed by a combination of rescalings of the variable £ 
and the energy. Using this freedom and starting from non-vanishing a and (3 one 
can always reduce them to "standard" a = (3 = —2, see below. The parameters 7 
and 5 remain free. 

Requiring the matrix (§) to have non- vanishing eigenvalues leads to a constraint 
that neither both parameters (0,7) nor both {(3,5) can be put to zero. One of 
the parameters in each pair can vanish, however. For instance, if a = the gen- 
eral elliptic potential degenerates into a polynomial potential. Thus, the example 
presented in Eq. ([]]) is nothing but a degenerate case of (§) - it corresponds to 
a = 0, 13 = 5 = -2, 7 = -(2j + l). 

Needless to say that j is an additional free parameter taking a discrete set of 
values. Thus, we deal with the three-parameter family of potentials: two continuous 
ones and one discrete. 

3. Elliptic potentials - special case. Prior to considering the general QES 
potentials with the ER symmetry we find it illuminating to discuss a few represen- 
tative examples. We start from 

a = (3 = -2 } 7 = -(8z/ + 6j + l), 5=-(2j + l), (10) 

where v is a constant. Since v is free, so is 7; the parametrization of 7 above, 
in terms of v and j, will be considered "standard". Physical arguments (e.g. the 
stability of the potential) require v to be non-negative. The parameter 5 is fixed for 
the time being. Later on we will let 5 vary too. 

The physical variable x in Eq. (|5]) is determined by the inversion of the equation 

lr The representation (|J) is a sufficient but not necessary condition. At some specific values of n 
there may exist QES systems with the ER symmetry which do not fall in the class of the systems 
we built. 
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(see |], |3|]). Equation ([□]) has solutions —V(x) and \jV{x) where V is the Weier- 
strass function. One could use either of them; the second solution is more convenient 
for our purposes. Thus, 



i[x) = —— g 2 = A, 53 = 0, (12) 
V[x) 

where #2,3 are the invariants of the Weierstrass function. For the time being it is 
assumed that £ £ [0, 1] and x £ [0, 2*] where 

rfS = V5T(5/4) m 
io 2v^F r(3/4) 1 ; 

Later on this constraint will be relaxed. Equation (0) maps the interval [0, a:*] onto 
[0, 1]. The function £(x) is double-periodic in the complex plane, with the periods 
2x* and 2ix t . Thus, under our choice of parameters, the parallelogram of periods of 
the Weierstrass function becomes square. The symmetry of the square immediately 
translates in the ER symmetry of the quantal problem at hand. We will use the fact 
that 

£(-*)= Z(ix) = -£(x) (14) 

stemming from the properties of the Weierstrass function with the above periods. 
The expansion of at x = runs in powers of x 2 . 

The phase a(x) of the gauge transformation and the "gauge potential" A(x) are 



1 f (cLA/dQ - 2B 
"(*)= Z J di 



-i/ln(l-£ 2 ) , (15) 



and 

1 {dA/dg) - 2B 4u^ 2 



A 



2 4^2A 



(16) 



As a result, we get the following potential in the Schrodinger operator 



V(x) = \ 4z/( f ] )e - (8f + 2j + I6uj + 6u)t) . (17) 

(The general formula for calculating the corresponding QES potential in the case at 
hand reduces to 



where C(£) is defined in Eq. (^j).) 

This Schrodinger problem is quasi-exactly solvable and can be considered beyond 
the original interval [0, x*]. For v = we deal with the periodic potential defined on 
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Figure 1: The periodic QES potential of Eq. ( |TTD at v — and j 

V(X) 60 r 




Figure 2: The QES potential of Eq. (|T7| ) at v = 1 and j = 2 defined on the interval 
— < x < x*. The potential has a double- well shape. 



the entire x axis (Fig. 1), which is akin to the Lame problem [10|. If v > 1/2 the 
potential is singular at x = ±x* (Fig. 2), and the problem is defined at x € (—a;*, x#). 
The condition v > 1/2 is necessary for stability. The Hamiltonian changes sign under 
the transformation x — > ix, 



H -y-H. 



x 



IX 



as follows from Eq. (fH]). The eigenfunctions {V'e, V'-b} m each pair interchange. 
The ER symmetry is explicit. Let us consider separately two cases, 
(i) v — 0. The potential and the eigenfunctions take the form 



V(x) 



2j(4j + l) 
V{x) 



, if>(x) = P 2j {£) 



where Pij is a polynomial of degree 2j. At positive E the spectrum is continuous, 
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while the counterpart wave functions ^-e with negative energy eigenvalues corre- 
spond to the boundaries of the Bloch zones; all these eigenf unctions are periodic. 
A very similar system, with a different coefficient in front of 1/V(x), emerges at 
v = 1/2. 

(ii) v > 1/2. The potential has a double- well form (Fig. 2). The singularity at 
x — ► is* is of the form (x±x*)~ 2 . The wave functions must vanish at x — ir*, the 
spectrum is discrete. The algebraic sector includes the ground state and 2j excited 
states symmetric under x — > —x. If j is integer, one level lies exactly at zero, j 
levels are below and above zero, respectively. If j is half-integer, (2j + l)/2 levels 
lie below zero and the same number above. 

4- Generic elliptic potentials with the ER symmetry. To proceed to the 

general case we invoke the only remaining freedom and let S in Eq. (H) float. The 
following parametrization will be used: 

<5 = -(2j + l)-/z, (19) 

while a, (3 and 7 are the same as in Eq. (|IUD. If /J 7^ Eqs. ( [T5D - fllTD are modified 
as follows: 

a(x) = LuHl-e)-^n T ^-\ , (20) 

A J±2L-J^\ , (21) 



and 



4e 



2^ 
1-e 



(22) 



where ^(x) is the same as before, see Eq. ([12]). By inspecting this potential one 
concludes on physical grounds that 

8u-A>/2>2. (23) 

Since the potential ( p2]) is singular at x — >■ (it explodes as 1/x 2 ) the problem is 
defined for \i > 2 on the interval (0, x*). If /i = 2 there is no singularity at x = 0, 
and the problem is defined on the interval (— £*,£*). The potential is depicted in 
Fig. 3. 

5. x — > ix, H — * —if. Since the Hamiltonian we built possesses this 

property one might ask why the ER symmetry is realized only in the algebraic 
sector rather than for the whole spectrum. From Fig. 2 it is quite obvious that 
the whole spectrum cannot have this symmetry - the states at positive E stretch 
indefinitely, while for negative E the lowest level is the ground state. Although the 
answer to this question is rather obvious, an explanatory remark is in order. 
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Figure 3: The QES potential of Eq. (|22"D at v = 1, // = 3 and j = 3 defined on the 
interval < x < x*. 



At arbitrary the second order differential equation Hip = Eip has two linearly 
independent solutions, let us call them ipi t 2- Generically both are non-normalizable. 
For E = E n one of the solutions (say, ipi) is normalizable, the other is not. Gener- 
ically, for arbitrary n, the transformation x — * ix connects a normalizable solution 
at positive E to a non-normalizable one at negative E. The latter does not lie in the 
physical Hilbert space, and there is no physical symmetry E «-> —E. However, if n 
belongs to the algebraic sector the transformation x — > ix does not generate a non- 
normalizable solution, since the phase factor a is invariant under x — > ix. Rather, 
in this CclSC X ^ IX connects a normalizable solution at positive E to another nor- 
malizable solution at negative E. Both belong to the physical Hilbert space, and 
E <-> — E is a valid symmetry. 

The argument above is somewhat simplified, we cut corners. Although the con- 
clusion is perfectly valid, the careful treatment would require explicit introduction 
of the Stokes lines and consideration of sectors in the complex plane. It is important 
that in the problems with the ER symmetry under consideration, one does not jump 
from one branch onto another in the process of the analytic continuation from the 
purely real to purely imaginary values of x. In this respect the situation is different 
from that discussed in |TI| , where the emphasis was on problems with the branch 
intertwining and the leaps from one branch onto another. 

6. The zero mode. If j is integer, there always exists a solution of the 

equation Hip = Eip with the vanishing energy eigenvalue, E = 0. The corresponding 
wave function ipo(£,) contains only even powers of £, so that it is invariant under 



x 



ix: 



it can be given in the closed form, 



MO = E 



J! 



k=0 



[7 - otj][l ~ Mj - 2)]..- [7 ~ Mj - 2fc + 2)] 

k\{j - k)\ [(j - l)/3 - 5)[(j - 3)/? - S]...[(j -2k + l)/3 - 5] 



e k . (24) 
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7. Quasiclassical vs. weak coupling expansions. Consider the lowest and 

the highest levels in the algebraic sector in the limit j > 1. We will denote them 
by i[)-Eo and i/jeq, respectively. The highest level is highly excited, and as such, 
is amenable to the quasiclassical treatment. The quasiclassical calculation of the 
energy of ipE was carried out f\ in Ref. M. On the other hand, ip-E corresponds to 
a system at the bottom of the well. This system is close to the harmonic oscillator, 
with a weak anharmonicity. One can develop the standard perturbation theory. The 
quasiclassical expansion and the weak coupling expansion have one and the same 
parameter, and coincide term by term, up to the overall sign |12| . 



Although the assertion above is quite general and refers to all QES problems with 
the ER symmetry, we will elucidate it using the simplest example. This will allow 
us to avoid bulky formulae. The generalization to the general case is transparent. 

As was noted in || , the simplest QES problem with the ER symmetry is that of 
Eq. (|I|). It is convenient to introduce a slightly different notation, 

H = -2T°T- - (2j + l)T~ - 2T+ , 



X 6 K 

Y~2 X 



.2 



■4>{x) = E^){x) , (25) 



Id* 
2~dx~ 2 

where 

K = 8j + 3 = 3,7,11,15, .... 

We are interested in the limit k — > oo. 

At large k the depth of the double-well potential becomes large and the well 
width small. The minima of V(x) lie at 

x = ±x , x = (k/3) 1/4 . (26) 



Near, say, the right minimum 

V(x) = -i^=- + 2k (x - x ) 2 + 10 i^-j (x - x ) 3 + ... , (27) 

where the ellipses denote quartic and higher order terms; similar expansion is valid 
near the left minimum. From Eq. fl2~7|) it is easy to find the ground state energy in 
the form of an expansion in 1/ k. Indeed, if we neglect exponential terms of the type 
exp(— Ck) where C is a positive constant, arising due to tunneling from one well into 
another 0, the ground state level can be considered as that of the harmonic oscillator 

2 Note that the unnumbered equation after Eq. (12) in ||] contains errors in signs. These errors 
miraculously combine with another error - the parameter of the quasiclassical quantization n was 
taken to be (k + l)/2 in Ref. J^], while actually it is (k — 3)/2 - to annihilate each other. The final 
expansion for the energy E presented in Eq. (13) of Ref. Q is perfectly valid. Note that our k 
corresponds to 4 J — 1 in § . 

3 The impact of tunneling and the issue of the analytic structure in k are discussed in more 
detail in Sec. 8. 
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slightly perturbed by cubic, quartic, etc. terms. The leading term in the ground 
state energy E is just the classical energy of the particle at rest in the minimum, 
i.e. at x , 

(kW 2 



3^/3 



The next-to-leading term is the zero-point oscillation energy of the harmonic oscil- 
lator, 



2 



A/2 



Then come the corrections due to the anharmonic terms in the potential (|2"7|). The 
first order correction due to the cubic term obviously vanishes. Therefore, the 
next term in the 1/k expansion of E comes from the quartic term in Eq. (|27|) 
(treated as a first order perturbation), plus the second order perturbation generated 
by 10 (k/3) 3/4 (x - x ) 3 : 

135 1 275 1 

~32~ 3V3K 1 / 2 ~ ~32~ 373^2 ' 

respectively. 

Assembling all these terms together we arrive at the following expansion for the 
ground state energy — Eq: 



-En = - 



3^ 35 



0(k 



(2* 



So far only the lowest level was discussed. What can be said about the highest 
level in the algebraic sector? 

The ER symmetry implies that the last level belonging to the algebraic sector 
has the energy 



En 



I) 



3/2 



1 



3^ 



35 ^ 



K 



in 2 



(29) 



Being considered as an excited state from the full set of states of the Hamiltonian 
(p5j), this level should have been labeled by (k — 3)/2. Indeed, i])-e is the ground 
state, then comes the first P-odd state, the first P-even excitation, etc. The ground 
state and 2j P-even excitations belong to the algebraic sector. The last (P-even) 
state from the algebraic sector has 



k — 3 



n 



(30) 



where n is the number of zeros in the corresponding wave function. 

Let us discuss now how the very same expansion for E emerges in the WKB 
approximation ||. It is instructive to start from the leading WKB approximation. 
Bohr and Sommerfeld's quantization rule at large n implies 



pdx 



nn 



K 



(31) 
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where 



P 



V2E - x 6 + kx 2 



and a and b are the turning points. We check that Eq. 
(ac/3) 3 / 2 . It is convenient to rescale the coordinate x, 



x = 2 l/2 3 -l/4 ft .l/4 2/> 



(32) 

lj) is satisfied at E = E = 

(33) 



Then 



p(E ) = V2 (£) ^ yfl - V + 3y 2 = V2 (£) v/(l-y 2 )(l + 2^ . (34) 



3/ 



3/4 



3/ 



At E = Eo the expression for E 1 — factorizes, and the integral / poke which in general 
is representable through elliptic functions in fact reduces to elementary functions || . 
Thanks to factorization we immediately see that the turning points are at y — ±1, 



j b V dx = -K j o dy(l + 2 y yi-y* = —, 



(35) 



q.e.d. 

The first correction in the quasiclassical expansion can be calculated as easily as 
the leading term. Indeed, at this level the only change to be done is the substitution 



1 K 

n — » n H — = 1 

2 2 



in the WKB quantization condition ([H]), and 

• \ 3/2 



E — > En 



(36) 



(37) 



where C\ is a numerical coefficient, to be determined from the quantization condition 



p(E )dx 



--If. 
2 



(38) 



Now p(En) takes the form 



p(E ) = V2 (if 4 



1 — — — % 6 + 3y 2 

K 



y/2 ('? 



\3/4 



l_ y2)(1 + 22/2) 2_^ 



2^ J(l-y*)(l + 2y*) 



+ 



(39) 



We have already checked that the O(k) term in Eq. ( j38| ) (it corresponds to keeping 
the first term in the square brackets) implies Eq = (k/3) 3//2 . Matching of the 0(k°) 
term in Eq. ( p5| ) (it corresponds to the second term in the square brackets) yields 



d = 3V3, 



(40) 
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in full accord with Eq. (|29|). 

Next-to- leading corrections in the quasiclassical expansion are calculated too 
||, see also footnote 1 above. The third and higher terms in the expansion require 
certain modifications of the WKB quantization condition which go beyond Eq. (j3£j). 
From what we already know about the QES systems under consideration, it is clear 
that the 1/k expansion of Eq obtained through WKB must match the weak coupling 
expansion. Six terms in the quasiclassical expansion of E were found in Ref. |§. 
Although it was expected, it was amusing to observe the coincidence with the first 
six terms in the weak coupling expansion. 

8. High-order behavior of the expansion. It goes without saying that 

the weak coupling expansion ( p8| ) is asymptotic. This is due to the possibility of 
the "leakage" from the right to the left well. The high-order terms are factorially 
divergent and of the same sign. The behavior of the high-order terms in the 1/k 
series for the ground state energy is determined [H| by the action of the instanton, 



the classical trajectory connecting the left and right minima in the Euclidean time. 
Let the instanton action be Sqk, where So is a number which we will calculate 
shortly. The Borel-resummed expression for the ground state energy has the form 

W vSw-" (4l) 

where the principle value prescription applies. The imaginary part of the integral is 
canceled by the imaginary part coming from the instanton-anti-instanton transition, 
which, in turn, is proportional to exp(— 2Sqk). The condition of cancellation 



fixes the denominator of the integrand. Expanding Eq. ( f4"T|) in 1/k we find the 
high-order tail of the Eq expansion, 

- Eq ~ -k'I 2 Y n\ - } - , (42) 

where no is an integer large enough for the asymptotics to set in. The instanton 
action is readily calculable, 



rxo 2k 3 / 2 /kX 1 ^ 
S k = J dx^ -^=- + x 6 - kx 2 , x = \^-J , (43) 

from where we obtain 

S = In i±^5 ~ 0.658 . (44) 

The ER symmetry and Eq. ( [12] ) imply that the very same factorial divergence 
is inherent to the quasiclassical expansion for energies of the highly excited states. 
Certainly, this phenomenon is known in the literature [ITJJ]. We find the argument 
above to be an illuminating way of demonstrating the asymptotic nature of the 
quasiclassical expansion. In fact, it is likely that the asymptotic regime starts quite 
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early. Indeed, the first five coefficients in the quasiclassical expansion can be inferred 
from Ref. ||. Denote the coefficients in front of l/(2SoK,) n by C n . Then, from Eq. 
(13) of this work we get 

C 3 ^6.57, C 4 ~20.2, C 5 «117, 

to be compared with the asymptotic prediction (|42|) 

C 3 = 3! = 6 , C 4 = 4! = 24 , C 5 = 5! = 120 . 

Barring the possibility of a coincidental proximity, we conclude that no can be as 
low as three. 

The parameter k is related to the cohomology parameter and is quantized. The 
nature of the 1/k expansions is closely related to the singularity structure in the 
complex k plane. In discussing this structure one should exercise caution, since 
the analytic continuation is performed from a discrete set of points, k = 3, 7, 11,.... 
This is one of the reasons why the singularity structure in the complex k plane turns 
out to be totally different from that discussed in earlier works JllJ, devoted to the 



analytic continuation in continuous parameters in the QES problems. There are also 
some other reasons responsible for the distinctions, e.g. k appears as a coefficient 
of a sub leading term in the potential, which is important. We do not dwell on this 
issue here, since it deserves a dedicated analysis. 

The quasiclassical quantization and the associated expansion imply k to be inte- 
ger (more exactly, k — 3, 7, 11, ...). At the same time, the weak coupling expansion 
(P5|) is the same independently of whether or not n e {3, 7, 11, ...}. It holds for any 
sufficiently large k. Both expansions coincide order by order, to any finite order; 
yet if k ^ {3, 7, 11, ...} the physical ER symmetry is absent, there is no reason for 
the coincidence of the absolute values of energy. This means that the factorially 
divergent weak coupling series and the quasiclassical expansion, presented in the 
square brackets in Eqs. (0), (0), respectively, define, generally speaking, two dis- 
tinct functions, despite the fact that the expansions per se are identical, order by 
order. The difference between these two functions is of the type sin(7r/t) exp(— Ck); 
it vanishes at k — 3, 7, 11, .... For these and only these values of k, making a full 2ir 
circle in the complex plane around k = 00, starting from a positive k and returning 
to the very same point, we smoothly interpolate between the lowest and the highest 
levels in the algebraic sector; their positions interchange. 

9. The ER symmetry in the finite difference problems. We have to mention 
that the energy reflection symmetry appears also in quantum-algebraic problems 
with the Hamiltonians built from finite difference operators (such problems naturally 
emerge in solid state physics). In order to display this property let us consider, for 
instance, the dilatation-invariant discrete operator defined as 

Ma = fJ %^f. («) 
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also known as the Jackson derivative (see e.g. |LoL |16f). Here q is a complex num- 
ber. In the limit q — > 1 the Jackson symbol obvioulsy goes into the conventional 
derivative. 



Now, one can easily introduce [17] a finite-difference analog of the algebra of 
the differential operators (Q) based on the operator D^, instead of the continuous 
derivative (for a discussion see J?|) 

T+ = {nR - , f° = -n + ^ , f~ = D ( , (46) 

where 



is the so-called q number, and 



l-q n 

\n\ = 

1 1 l-q 

{n}{n + 1} 
{2n + 2} 



It is easy to check that the operators (^) obey the commutation relations of 
the quantum algebra s^q for any value of the parameter j = n/2 (see e.g. [JIB] ). If 
j is a non-negative integer, the finite-dimensional representation (§) of the algebra 
(|46|) exists; it is irreducible when q is not a prime root of unity. The same line of 
reasoning which we followed to demonstrate the ER symmetry of the Hamiltonian 



(§) can be used in the case of the finite difference generators (f|6|). In this way we 
arrive at the conclusion that the discrete Hamiltonian 

H = af + f° + (3f°f- + -ff + + 5f- = 

A(ODj + B(OD^ + (7(0 , (47) 

possesses the ER symmetry. Here a, (3, 7 and S are numerical constants, and A, B, C 
are polynomials of the third, second and first degree in £, respectively, 

A{0 = + K , B(£) = [a({n} + n - 1) - 7 ]£ 2 + {6- P{n}) , 

£(£) = {n}( 7 -an)£, (48) 

where n = 2j. 

It is remarkable that a particular form of this quantum-algebraic Hamiltonian 
(with a slightly different definition of the discrete derivative) appears in the Azbel- 
Hofstadter problem of the electron motion on the two-dimensional lattice in the 
transverse constant magnetic field |T9|, EDI. In this case the parameter q is a prime 



root of unity; it is related to the magnetic flux through the lattice plaquette (the 
flux is given by a rational number with an even denominator). 

10. Comment on the literature. In Ref. [^TJ a certain "duality" trans- 

formation was suggested for the QES systems which inverts the signs of all levels 
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belonging to the algebraic sector and, simultaneously, changes the form of the poten- 
tial in a concerted way. It was observed that the potential (P is self-dual. Thus, the 
ER symmetry of the Schrodinger problem ([TJ) was rediscovered. It was noted then 
that the quasiclassical treatment of the QES problems should be qualitatively differ- 
ent from that of "conventional" problems, where there is no (quasi)exact solvability. 
The corresponding remark in |21| is rather vague, and we feel that an explanatory 
remark is in order here. 

Suppose the wave functions of a quantal system are treated in the WKB approx- 
imation. The WKB asymptotics, being considered in the complex x plane, contains 
singularities at the points where the classical momentum vanishes. The Stokes lines 
are attached to these points; they divide the complex x plane into several sectors. 
The appropriate WKB expression for the wave function in the given sector, when 
analytically continued across a Stokes line, may or may not match the appropriate 
WKB expression in another sector. In other words, distinct asymptotics may apply 
in the different sectors in the complex plane. This is a general situation. In the 
QES problems, for those levels that are determined algebraically, the wave function 
is analytic everywhere except infinity []. One and the same asymptotics remains 
valid in all sectors; one can freely do analytic continuations across the Stokes lines. 
The singularities of separate parts of the WKB expressions for the wave functions 
are superficial; they cancel when all parts are assembled together. This property is 
well-known in the harmonic oscillator, it extends to all QES systems, however. 

The observation above belongs to A. Vainshtein. He pointed out that the re- 
quirement of cancellation of these apparent singularities can be used in order to 
generate QES potentials. This requirement acts as a substitute of the algebraic 
structure within the Lie-algebraic approach. 
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